Towards Determining the energy of the UHECRs observed 

by the ANITA detector 

Konstantin Belov for the ANITA collaboration 

UCLA Department of Physics and Astronomy / 154705, 475 Portola Plaza, Los Angeles, CA 90095, USA 

Abstract. The Antarctic Impulsive Transient Antenna (ANITA) is a balloon-borne radio experiment designed to discover 
ultra-high energy cosmic neutrinos. The ANITA detector has completed one prototype and two full-scale flights above the 
Antarctic continent. Two direct and fourteen reflected cosmic ray events of ultra-high energy were observed during the first full 
scale flight and several others in the second flight. We present a Monte Carlo technique and analysis developed to determine 
the energy of the primary cosmic ray particles from the ANITA data. 
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1. INTRODUCTION 

Cosmic ray particles are constantly bombarding the 
earth's atmosphere. The energy of some of the cosmic 
rays can exceed 10^^ eV [1], making them the carriers 
of valuable information about the Universe and particle 
interactions at ultra-high energies. These ultra-high en- 
ergy cosmic rays (UHECRs) have been the object of in- 
tensive study since the 1938 discovery by Pierre Auger 
of extensive air showers, cascades of secondary particles 
in the atmosphere caused by UHECRs. Despite signifi- 
cant progress in research, the origin, acceleration mech- 
anism and chemical composition of UHECRs are still 
unknown. The flux of the cosmic rays falls rapidly with 
energy, making direct observations of UHECRs impossi- 
ble. 

Indirect observations employ the earth's atmosphere 
as a giant calorimeter to observe secondary particles on 
the ground using muon counters, or to observe air fluo- 
rescence or Cherenkov light produced by the extensive 
air showers using special UV light sensitive cameras. 
The latter observations can only be conducted during 
moonless nights, reducing the duty cycle of the detec- 
tor to less than 10%. A hybrid technique combines air 
fluorescence detectors and ground counters, improving 
data statistics as well as arrival direction resolution, pri- 
mary particle energy and the depth of the shower maxi- 
mum (Xfnax) reconstruction. This latter parameter of air 
shower development in the atmosphere can be used to 
obtain chemical composition and particle cross-sections 
at ultra-high energies [2, 3, 4]. The ground array and air 
fluorescence techniques have advanced significantly over 
the last several decades, but radio detection still has sev- 
eral significant advantages, including lower deployment 
and operational costs combined with a 100% duty cycle. 
Unfortunately, the interest in radio detection diminished 



after the initial experiments in the 1960's [5]. The lack of 
understanding of the radio production mechanism at that 
time resulted in somewhat limited success. 

A new breed of ground radio detectors has appeared 
in the last several years [6, 7, 8, 9] complemented 
by the Antarctic Impulsive Transient Antenna (ANITA) 
balloon-borne detector, which reported the first observa- 
tion of UHECRs in the 200-1200 MHz frequency band 
in a far zone [10]. The ANITA instrument is a self- 
triggering detector, which puts it apart from ground radio 
arrays usually located near an existing cosmic ray tele- 
scopes that supply trigger information. With the devel- 
opment of Monte Carlo (MC) simulations incorporating 
precise physics [11, 12], it becomes possible to obtain 
the energy of the primary cosmic particle as well as in- 
formation about the air shower development from radio 
data. This makes radio observations more attractive as a 
complement to the existing cosmic ray telescopes or as a 
stand alone independent observation technique. 

A special approach, however, is needed to determine 
the energy of the primary particle and the air shower 
development parameters from the radio data. We present 
a technique to obtain the energy of the UHECRs using 
MC simulations in conjunction with the total power and 
spectral density of the observed radio pulses emitted by 
the extensive air showers. 

The radio emission from an extensive air shower in- 
creases with primary energy. Due to relativistic ampli- 
fication, the maximal power will be emitted along the 
Cherenkov angle of the shower. The proposed technique 
employs the spectral power density to identify the shower 
core position with respect to the receiver. To obtain the 
primary particle energy, one has to combine it with the 
total broadband power measured by the receiver. Because 
this technique relies only on measuring the radio emis- 
sion, it proves to be a powerful method for measuring 



UHECR energies not only with balloon-borne or space- 
borne receivers, but also with ground-based detectors. 
With the ability to reconstruct the Cherenkov ring, the 
ground-based detectors can also measure the depth of the 
shower maximum, X^ax, opening a possibility to study 
the cosmic ray chemical composition and particle cross- 
sections at ultra-high energies using the radio data. 



2. ANITA DATA 

ANITA is a radio detector designed to observe ultra-high 
energy cosmic neutrinos [13], see Figure 1. An ultra-high 
energy neutrino interaction in ice will emit a radio pulse 
due to the charge asymmetry in the cascade of secondary 
particles, the so called Askaryan effect [14, 15]. In case 
of a neutrino skimming the ice, the radio pulse can be de- 
tected from above. Flying some 37 km above the Antarc- 
tic ice, the ANITA detector is sensitive enough to pick 
up weak radio pulses from hundreds of kilometers away 
and to reconstruct the event location to exclude events 
with anthropogenic origin. The high altitude provides a 
large detector aperture while the deep Antarctic ice, with 
about 2.5 km average thickness, provides a large detector 
volume. 

The ANITA prototype flew in the 2003-04 Antarctic 
season, followed by the ANITA I flight in 2006-07 and 
the ANITA II flight in 2008-09. The third ANITA flight 
is scheduled for the 2013-14 Antarctic season. The two 
ANITA flights did not find ultra-high energy neutrinos 
above the expected background, and a limit on the neu- 




trino flux is published [16]. The neutrino signal is ex- 
pected to have a vertical polarization due to the event 
geometry. However, 16 isolated horizontally polarized 
events were identified in ANITA I data, see Figure 2. 
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FIGURE 2. Map of locations of detected reflected (red di- 
amond) and direct (black square) UHECR events superim- 
posed on a microwave backscatter amplitude map of Antarctica 
(Radarsat) within ANITA field of view (dash-dotted line) [10]. 
For the direct events the closest approach to the surface of the 
RF signals is shown. The scale on the right is relative amplitude 
of the backscatter signal. 

These events cannot be associated with ultra-high en- 
ergy neutrino interactions, and are caused by UHECRs. 
The 14 events reconstructed on the ice have the same 
electric field polarity, see Figure 3, while the two events 
reconstructed above the horizon have opposite polarity, 
which is consistent with the geosynchrotron mechanism 
of radio emission. The radio signal from the 14 exten- 
sive air showers, beamed up front into a cone, is reflected 
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FIGURE 1. ANITA I before the flight. Photo by Jeff Kowal- 
ski. 



FIGURE 3. Waveforms of 14 ANITA horizontally polarized 
reflected events, shifted and plotted on the same scale. The 
instrument response is deconvoluted. Plot by S. Hoover. 
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FIGURE 4. ANITA UHECR detector concept. 



from the ice surface and the polarity of the electric field is 
flipped. An additional argument supporting the geosyn- 
chrotron mechanism as cause for the radio signal is that 
while the ANITA detector did not trigger on the weak 
vertically polarized content of these events, the residual 
electric field in the vertical polarization of each event 
correlates with that expected at the reconstructed event 
location [10]. In fact, ANITA is a sensitive UHECR de- 
tector with a very large aperture, see Figure 4. To ob- 
tain the energy of the UHECRs from the ANITA data we 
have to use MC simulations, taking into account radio 
frequency (RF) signal propagation in the air, reflection 
from the rough ice surface and instrument response. 

3. MONTE CARLO SIMULATIONS 

CORSIKA code version 6.960 [17] was used for the 
extensive air shower simulations for this work, and the 
RF emission from the air shower was calculated using the 
CoREAS plug-in for CORSIKA [18]. The simulations 
were performed for proton primary particles at different 
energies from 10^^ to 10^^ eV using the actual event 
geometries reconstructed by the ANITA detector. The 
atmospheric parameters were chosen for the South Pole 
January atmosphere, and the magnetic field at the event 
locations were used. The MC simulations also take into 
account the effects of the Fresnel reflection, roughness, 
antenna response and the RF propagation to the payload 
as well as the realistic index of refraction, changing with 
the atmospheric depth. The electric field intensity on the 
surface of the ice for one of the simulated ANITA events 
is shown on Figure 5. As expected for the real data 
in case of the index of refraction greater than one, the 
simulation shows that the radio emission forms a ring, 
actually an ellipse, on the ice surface, corresponding 
to the Cherenkov cone projection on the ground. For 
observers inside the ring, the time is reversed, that is, the 
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FIGURE 5. Horizontal polarization component of the elec- 
tric field on the ice surface at the time of maximum electric 
field magnitude for a simulated ANITA event. Coordinates are 
relative to the shower core. 

observer sees the end of the shower first. For observers 
located somewhere on the ring, the shower develops 
instantaneously and the emitted signal is amplified due 
to relativistic effects. The shower core is in the focus of 
the ring, and the ring geometry depends on the air shower 
zenith angle and Xf^ax- 

4. ENERGY OF THE UHECRS 

The ANITA detector sees the signal reflected from a 
location on the Cherenkov ring. Because the detector 
does not see the signal from the whole ring, the total 
power of the received radio pulse alone cannot be used to 
determine the energy of the primary cosmic ray particle. 
This is true not only for a balloon-borne detector like 
ANITA, but also for a ground array of radio antennas. 
Indeed, an observer on the Cherenkov ring can see a 
very strong signal while an observer inside or outside 
of the ring will only register noise for the same cosmic 
ray event. Figure 6 shows the signal intensity received 
by the payload as a function of primary particle energy 
and radial distance from the shower core for a series 
of simulated events of different primary energy. The 
geometry of one of the 14 ANITA cosmic ray events was 
chosen for this simulation. 

Using only the total power of the received radio signal 
to determine the energy results in a large uncertainty, be- 
cause the received signal may have been reflected from 




FIGURE 6. Total power in 200-1200 MHz band in units 
of dBm received by the detector for a series of simulated 
events of different primary energy. An actual ANITA event 
geometry was chosen for this simulation. Signal propagation in 
the atmosphere, ice reflection effects and instrument response 
are accounted for. 



any location on the ring. While there is an obvious "az- 
imuthal" coordinate degeneracy of the reflection point 
seen by the ANITA detector, it is possible to measure the 
radial distance from the shower core. The signal fluctua- 
tions along the Cherenkov ring are significantly less than 
those in the radial direction and can be included into the 
systematic error. 

In addition to finding the total power, a broadband 
measurement like that made by ANITA can also deter- 
mine the power spectrum of the received radio impulse. 
The spectrum can be used to determine the reflection 
point radial distance from the shower core, the x-axis co- 
ordinate on the plot in Figure 6, thus limiting the uncer- 
tainty on the primary particle energy. Figure 7 shows the 
power spectrum density for the same simulated event at 
different distances from the shower core, together with 
the actual event data recorded by the ANITA detector. 
Note the 5 m scale used for the simulation. 

The power spectrum for the simulated data and 
ANITA events is quantified using a spectral ratio, R = 
P2/P1, where Pi is the total power in the 350-600 MHz 
and P2 is the total power in the 650-850 MHz frequency 
range, respectively. It is also possible to use the spec- 
tral index for this purpose, but for noisy radio signals the 
spectral ratio is usually more stable. The MC simulations 
have shown that, as expected, the spectral ratio R does 
not depend on the primary particle energy, but changes 
significantly with the distance from the shower core, 
reaching its maximum in the middle of the Cherenkov 
ring, see Figure 8. The maximum of the spectral ratio 
corresponds to the hardest radio spectrum observed from 



FIGURE 7. Power spectrum density of the radio emission 
from the simulated air shower and the ANITA event (black 
line). The smoothed ANITA data (green line) is also shown. 



the air shower. On the Cherenkov ring the signal appears 
coherent at all frequencies, but the coherence diminishes 
faster for higher frequencies as the observer moves away 
from the middle of the Cherenkov ring, leading to a 
softening of the spectrum. The radial distance from the 
shower core is fixed based on the spectral ratio analysis, 
which significantly reduces the primary particle energy 
uncertainty, see Figure 6. The total power of the ANITA 
event, of which the geometry was used for this simula- 
tion, is -85.69 dBm. Combining the total power of the ra- 
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FIGURE 8. Spectral ratio R for the ANITA event (green line) 
and the simulated event (black line) at different distances from 
the shower core. There is no signal far from the Cherenkov 
ring (0-220 m and more than 370 m from the shower core). 
Those regions are excluded from the analysis. Because of the 
noise floor at about -120 dBm/MHz, 20% of the total power is 
subtracted from the signal at the higher frequency band for the 
real data. The MC data does not include noise. 



dio signal, received by the payload, with the limits on the 
distance from the shower core, obtained from the spectral 
analysis of the radio signal, one can determine the energy 
of the primary cosmic ray particle. The uncertainty on the 
energy measurement depends on several factors includ- 
ing the uncertainties on the signal reflection from the ice, 
spectrum measurements and the systematic error on the 
radio signal measurement. 

It should be noted that it is not surprising that the 
ANITA detector sees the signal reflected at a distance 
from the focus that is centrally located within the 
Cherenkov ring. The detector was optimized for ultra- 
high energy neutrino observations, with the trigger sen- 
sitive to very short radio pulses that have broadband and 
relatively flat spectra. This trigger bias makes the detec- 
tor sensitive to the middle region of the Cherenkov ring 
only, and pushes the energy of the observed UHECRs 
lower by reducing the detector trigger aperture. 

5. DISCUSSION 

The radio spectra from extensive air showers produced 
by UHECRs provide the additional information needed 
to measure the energies of cosmic rays. For a balloon- 
borne or a space-borne radio detector the spectrum re- 
duces the energy measurement uncertainty, putting radio 
observations on a par with the air fluorescence and the 
ground counter techniques. 

The technique presented here can also be employed 
by ground arrays. For a ground radio array the spectral 
measurements allow for the full reconstruction of the 
Cherenkov ring using a smaller number of surface ra- 
dio antennas than those needed in absence of the spec- 
tral measurements. Indeed, the width, size and overall 
geometry of the Cherenkov ring on the ground varies 
widely, depending on the frequency range as well as on 
the extensive air shower zenith angle and the depth of 
the shower maximum. The dependence on Xynax opens an 
opportunity for studies of the cosmic-ray chemical com- 
position and strong interaction particle cross-sections. 
However, a large array of densely spaced surface anten- 
nas is needed to reconstruct the Cherenkov ring, while 
still covering an acceptable range of air shower param- 
eters with meaningful statistics. The spectral measure- 
ments provide the information about the location of the 
observer relative to the shower core, effectively reduc- 
ing the antenna density requirements, and, thus, the total 
number of antennas needed and the cost of the array. This 
opens a possibility of building radio arrays that comple- 
ment the existing detection techniques, or even operate 
as stand alone UHECR observatories. 
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